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The effect of Coulomb interaction upon superconductive proximity effect in disordered metals 
is studied, employing newly developed Keldysh functional approach. We have calculated 
subgap Andreev conductance between superconductor and 2D dirty film, as well as Josephson 
coupling via such a film. Both two qualitatively different Coulomb effects - suppression of 
the tunneling density of states and disorder-enchanced repulsion in the Cooper channel - are 
shown to be important at sufficiently low temperatures, In^ Tt > l/g, where g = Tia/e^ is the 
dimesionless sheet conductance of the film. 



1 Introduction 

Electron transport in hybrid superconductive-normal (S-N) systems at low temperatures is gov- 
erned by the Andreev reflection. Both flnite-voltage conductance Ga between superconductive 
and normal electrodes, and Josephson critical current Ic between two superconductive banks, 
separated by normal region, are determined by the Cooper pair propagation in the normal 
metal. The theory of Andreev conductance (without Coulomb effects) was developed in e. g.0'§i, 
whereas Josephson coupling was calculated (with the account of short-range electron interaction) 
ini. When normal conducting region is made of a dirty metal film, or two-dimensional electron 
gas with low density of electrons. Coulomb interaction in the normal region may lead to strong 
quantum fiuctuations which suppress both Andreev conductance and Josephson proximity effect. 

There are two different quantum effects which play major role in dirty superconductive 
systems. First of all, critical temperature of a uniformly disordered superconductive film 
is suppresed by disorder, and eventually vanishes Bi'i at the critical value of conductance 
gc = {2it)~^ In^ T^'i ^^^re r is the elastic scattering time and Tco is the BCS transition temper- 
ature. In this range of parameters weak-localization iH and interaction-induced corrections 
to conductivity of 2D metal are of the relative order of g^^ In <^ 1. This fact brings about 
the issue of a quantum Jj.e. T = 0) superconductor-metal transition (cf. as opposed 

to the usually assumed Hj direct S-I transition. The second effect that is important for hybrid 
superconductor-normal systems with tunnel barriers at the S/N interface, is the effect of tun- 
neling conductance suppression J2' zero-bias anomaly") due to slow relaxation of an extra charge 
after the tunneling event t9't3i§ll3; as noted intj, this is physically the same effect as Coulomb 
blockade of tunneling into a finite system. In a 2D system with a long-range Coulomb interaction 
this effect scales as In^ Whereas the first discussed effect is insensitive to the form of 
long-range Coulomb asymptotics, the zero-bias anomaly does crucially depend on it^and may 
be suppressed by screening of the Coulomb interaction by nearby external electrodes t3. 

In the present paper we discuss behaviour of the subgap conductance between superconduc- 
tor and thin dirty metal film, as well as Josephson coupling between two superconductors via 
such a film, assuming that relevant energy scale huj = max(ey, T) is such that < Inwr <^ g. 
The second inequality allows us to neglect weak localization and interaction-induced corrections 
to the film conductance, whereas the first one calls for a nonperturbative treatment of the zero- 
bias anomaly (ZBA) and Finkelstein effectsi'i0. The first of these effects was recently considered 
in £3, where renormalization of the tunneling DoS on the N side of N/S sandwich was calculated 
within replica functional method (cf. i). We prefer to avoid the useof replicas; instead, we 
have developed the functional integral approach based on the Keldyshlla representaion of Green 
functions for dirty super conductor scJ. An advantage of this method is that it allows to calculate 
nonequlibrium quantities, and does not contain any analitic continuation procedures. Recently 



this method was developed for normal metals incj; in many respects we follow that approach. 
Here W£ just present our main ideas and results, whereas full details can be found in forecoming 
paper cil. 



2 Keldysh cr-model 

In the simplest case of spin-independent interactions, and neglecting relativistic effects, we can 
write down an effective low-energy action for a dirty film of superconductor as EH 



S = —Tr 
4 



Diyqf + Ai{iT^dt + * + A)q\ + 2uTT$^a^$ + Tr (^o-^^Fo" V+ Tr A+cr^A 

J A 



Here Q = Qt^t' is a 4 x 4 matrix in the K^N space, which depends upon two time arguments t, t' . 
Pauli matrices in the Keldysh (K) 2x2 space are denoted by (yx,y,z^ where T^^y^z stand for the 
Nambu (A^) space. Matrix Q is subject to the usual a- model constraint = \. Operation Tr 
means taking trace over K®N matrix spaces, as well as over time and real spaces, v is the bare 
DoS per one spin projection at the Fermi level, D is the diffusion coefficient, (j) is the (fluctuating) 
electric potential, A is the BCS order parameter field, A is the interaction constant in the Cooper 
channel (its negative sign corresponds to attraction). Vo((7) = 27re^/g is the Fourrier component 
of Coulomb interaction in 2D. We use both vector and tensor (boldfaced) notations for bosonic 
fields in the Keldysh space, e. g. A = (Ai, A2)^ and A = AiCTq + A2(Ti, where Ai^2 are certain 
linear combinations of the A field on the forward and backward branches of the Keldysh time 
integration path. The field = (p — dfK, where K is the phase of the gauge transformation 
similar to the one proposed inE2l; in the superconductive case the original matrix field Q changes 
under this gauge transformation according to Qt,t' = e^^^^^^^Qt^t'e"*^*-* . Covariant space 
derivative VX = drX + {[TzdrK., X]. The last term in (||), denoted below as S^, describes 
electron tunneling across the barrier (cf. 0). Here Qs and Kg refer to the S side of the interface 
boundary F; the notation Trp means that the space integral is taken over the interface surface, 
7 is the (dimensionless) normal-state tunneling conductance per unit area of the boundary. 
Variation of the total action with respect to Q and Qs with the constraints = = 1 
leads (at = 0) to the dynamiclll Usadel equationil, together with the standard §1 boundary 
conditions. In the absence of the boundary term, the equilibrium saddle-point solution of the 
Usadel equation is given by Q = u • diag((5^, Q^) ■ u, where Q^^^") are retarded (advanced) 
Green functions, u = = o'z-'r cr^F, and F = f + /it^ is the generalized fermion distribution 
function, 0"+ = (ax + iay)/2. This representation suggests the use of a new variable, Q = uQu. 

Below we consider low-energy limit e <C |A|, so the superconductive matrix reduces to 
purely phase rotations: Qs = — iA/|A|. We will also assume weak tunnlling across S/N in- 
terface, which makes it possible to start (on the normal side of interface) from purely metal- 
lic saddle-point Qq = GzTz- Slow rotations of this trivial solution can be parametrized as 
Q = I'^GzTzS^ 1"^ by the matrix W subject to the constraint {W, ct^Tx} = 0, which is re- 
solved as W = I ^^^^ + ^f^-^ \ . The diagonal (off-diagonal) in the Nambu 

\^ Wo + WzTz WxTx + VJyTy j ^ 

space excitations, Wi and Wi with i = 0, z {i = x,y), correspond to diffusion (Cooper) modes in 
the metal. FollowingE^I, we choose K to be a linear functional of ip and require the vanishing of 
the term bilinear in W and $ in the cj- model action (|^). 

The term in the action, which is responsible for Andreev subgap conduction, comes from the 
averaging of over Cooperon modes: Sa = %C!a T^t^{Qs(^zTz)'^ , where Ga is the (dimensionless) 
Andreev conductance. At low energies {uj,T,eV) <^ Eth = fiD/L"^, L being the relevant length 
scale of the N region, and in the absense of quantum corrections, the known resulti = Gj^Rd 



is recovered, where is the total resistance of the diffusive metal region and Gt = 'J A is the 
total normal-state tunneling conductance, A being the area of the junction. In the geometry of 
rectangular N/S contact (normal metal film with the length Ly along the boundary with SC, 

(0) T 

and the length between SC and normal reserviour), G'X' = ^t;- At higher w and/or T, the 
length Lx should be replaced by Lefj(u;,T), which is equal to Leff(u;) = yj2D/\oj\ at huj ^ T, 
and to Lefr(T) = 0.95\/ D/2T in the opposite limit. 



3 Renormalization group 



Quantum fluctuations lead to renormalization of the Cooper-channel interaction constant A in 
the N metal and of the barrier transparancy 7; note that the sheet conductance g = 2vD is 
constant within our approximation. Renormalization of A comes about in the second order over 
the interaction term S^^q = — vri/Tr [n^uQ], which should be averaged over diffusion/Cooperon 
modes of the Q matrix field and electric potential fluctuations (j). This calculation can be done 
in the standard gauge with K = Q (cf. Ell for detailed discussion), the result coincides with the 



one known fromt 



1 ,2 . . 1 



A^; C = ln— , (2) 



dC, ^Ti'^g ' OJT 



where u is the running low- frequency cutoff of RG procedure, which stops eventually at u; = 
max(r, eV, E^h)- The last term in (|2|) is the usual BCS ladder contribution, which can also be con- 
sidered within RG approach, as coming from the second order term over Sa,q = — vrz^ Tr [uAnQ]. 
Eq. (^) has a locally stable infrared fixed point A^ = Xj^-K^. Integration of Eq. ^ with initial 

condition Aq < leads to Finkelstein's result for the supression: TcT = ( 1^^^ in(i/r'"or) ) ^ ^ ■ 
There are two sources of renormalization corrections to 7. The first one is due to Cooper- 
channel interaction; to find it one should average over Cooperon modes the product S^Sx, 
where Sx = (l/4)7r^z^A / dtcPrtr ax[Qtf — (tzQu)"^] is the result of Gaussian integration over A 
field in the action (||). The second contribution to 7 is due to fluctuations of the Coulomb- 
induced phase K(t), describing the ZBA effect; here we neglect similar fluctuations of the phase 
K-sit), assuming that voltage at the superconductor sideis fixed by external curcuit, and SC 
diffusion constant Ds ^ D (for the opposite situation cfBJ). Mean-squared fluctuations of the 
field K(t) can be expressed via the function p{uj) which generalizes a notion of "environmental 
impedance" introduced inB3 within phenomenological approach to the Coulomb suppression of 
Andreev conductance. We calculate p{u!) microscopically for all considered geometries, see 
below. The resulting RG equation for 7 reads 

At high frequencies uj > A one has instead dj/dd^ = —7 p/vr, as it should be for the tunneling in 
the absence of superconductivityi'0@. Eqs. (||), are sufficient to study low-energy behaviour 
of the Andreev conductance in a rectangular geometry. Another interesting problem is the 
subgap conductance between small SC island of size d and large metal film (cf. EJ), when the 
"diffusive" resistance Rd logarithmically depends on the relevant space scale. As a result, the 
Andreev conductance becomes itself subject to the RG transformation: 

^=A^^-^pGa at C>Ca, (4) 

where we assume, for simplicity, Ca ~ Cd = In 



4 Andreev conductance 



We start from the case of rectangular geometry of the contact, when it suffices to find renor- 
mahzed value ■^r by integration of Eqs. @, @, so that Ga = G'^^(7ij/7o)^. In the practically 
interesting (at < A) limit ^ L ^ Lgfj (a;) we find p{uj) = ^ In ^ ^^^^ (we assumed here 
L = Lx ^ Ly). The result for the Andreev conductance reads 

G|,L^ ,.e.pf-^,n^.l„^^ (5) 

9 i» [1 + A„/A, + (1-A„/A5)(ij/Ap.]' ^'g ^ El)'' 

where A„ is the Cooper-channel replulsion constant at the energy scale hT~^, Gta is the normal- 
state tunneling conductance at eV A. The last multiplier in the r.h.s. of Eq. (|5|) describes 
the ZBA effect, with the doubled coefficient in front of the log accounting for the fact that 
in this geometry extra charge spreads (after tunneling) over half-plane. Eq. (^) is valid at 
^ ^ Eth = In the opposite limit the result (^) can be used after replacements a; — > 

Next we consider geometry of small SC island siting in the middle of the N thin film of 
characteristic size L. Now we need to integrate the whole set of Eqs. (|2|)-(Q). We assume, for 
simplicity, that the island size d ~ .^a = yD/A. Thus at the energy scale uj <^ A extra charge 

1 In ^ 

2ng Leff(i^)' 



speads symmetrically over 2D plane and p{io) = In r ) ■ The result reads (assuming again 



^'^^^ = (A, + A„) + (A,-A„)(./A)-. • ^^H"^ - ' ^ J 

At lowest energies uj should be replaced by Et^ in Eq. (^). The whole effect of quantum fluctu- 
ations in this geometry is that initial " semiclassical" logarithmic growth of Ga with lo decrease, 
Ga oc In crosses over to the decrease by the "log-normal" law. 

Note that in both cases (|5|), (P) the ZBA effect is represented as a separate multiplicative fac- 
tor; the reason is that fluctuations responsible for the ZBA have frequencies uj ^ Dq^, and thus 
they are not mixed with the low-u; fluctuations respon sible for the A renormalization. Detailed 



dependences of the ZBA factor on V and T can found using formulae fromE3 with p{uj) em- 
ployed as the effective impedance. The power-law factors due to Finkelstein's renormalizations 
are determined always by the largest of the scales eV,T,Eth. 

5 Josephson proximity coupling 

The term in the effective action, which is responsible for the Josepshon proximity coupling, 
can be written in the form Sj = ^EjTrlQ^'^Q^pax], where superscrits and refer to two 
superconductive banks or islands. Using low-energy representation for Qs, and neglecting phase 
factors exp(iK(t)r^), one finds that calculation of the Josephson current with the use of Sj 
produces standard expression Ij = ^Ej sin{9i — 62) ■ The contribution of K(t) fluctuations 
factorizes as above, and is taken care of by the same multilicative factor as in Eqs. (|5|), (P), with 
UJ replaced by an appropriate inverse diffusion time. We start from the calculation of Ej in the 
absense of the ZBA. For both geometries discussed above, Ej can be expressed via the same 
function defined in the Fourier space, the zero-frequency Cooperon amplitude J[q)^ which obeys 
the RG equation ^'^^^''^ = •^'^-^-^r'l where (q = \n{A/ Dq'^). First we consider an example of two 
small SC islands of radius d, separated by the distance d. In this case Ej, which is given 
by the inverse Fourrier transform of J^{q), can be represented as Ej[R) = -^^^ ^'^^''^ ^ 1 ' 



Using renormalized 7(Cg) according to Eqs.(^J^, and adding the ZBA factor from (^) estimated 
at u = D/i?^, we find (with the total size L of the film being in the range R <^ L <^ Attuc'^R'^ ): 



where /3 = ^^q:^- The increase of the power- law exponent up to xj = 2 + 4Ag is due to 
Finkelstein's corrections, whereas additional log-normal decay factor is due to the ZBA. 

In the case of rectangulareeometry the discrete nature of diffusion modes in the N region 
should be taken into accounted while calculating Ej. The result (for ^ Ly ) is: 
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